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ABSTRACT : 

Two  algorithms  have  been  formulated  for  scheduling   n  jobs 
through  a  single  facility  to  minimize  the  number  of  late  jobs 
when  set-up  times  are  sequence  dependent.   The  first  is  a 
simple  matrix  algorithm  which  solves  the  problem  when  jobs 
must  be  processes  in  first-come,  first-served  (FCFS)  order. 
The  second  is  a  branch  and  bound  technique  which  arrives  at 
an  optimal  solution  with  no  restrictions  on  the  sequence  used. 
Both  algorithms  are  demonstrated  by  examples. 
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I.  INTRODUCTION 

The  problem  treated  is  the  n  job,  one  machine  scheduling 
problem  in  which  set-up  times  are  sequence  dependent.   Algorithms 
are  presented  which  minimize  the  number  of  late  jobs  for  two  var- 
iations of  the  problem. 

In  the  first,  jobs  are  restricted  to  be  processed  in  first- 
come,  first-served  (FCFS)  order.   After  presentation  of  this  al- 
gorithm it  is  applied  to  an  example  problem  and  then  proved  to 
produce  an  optimal  schedule.   The  following  section  draws  a  com- 
parison between  this  algorithm  and  Moore's  algorithm  which  produces 
optimal  schedules  when  set-up  times  are  sequence  independent. 

The  second  algorithm  produces  optimal  schedules  when  jobs 
may  be  processed  in  any  order.   A  branch  and  bound  technique  is 
developed  and  illustrated  with  an  example  problem.   Finally,  the 
amount  of  computation  required  by  this  algorithm  is  compared  to 
that  required  for  the  algorithm  of  Little,  et.al.,  to  solve  an 
n  city  traveling  salesman  problem. 

II.  PROBLEM  DESCRIPTION 

A  single  machine  is  used  to  process  a  known  set  of  jobs. 
Each  job  has  only  one  operation  and  its  processing  time  and  due- 
date  are  known.   The  machine  set-up  time  for  each  job  is  known 
but  is  dependent  upon  the  job  which  precedes  it.   A  job  once  begun 
is  processed  until  completion.   The  measure  of  performance  is  to 
minimize  the  number  of  late  jobs  when 


1)  Those  jobs  which  are  included  in  the  schedule  must 
be  processed  according  to  first-come,  first-served 
(FCFS)  discipline, 

2)  jobs  may  be  processed  in  any  order. 

The  nature  and  difficulties  associated  with  this  problem  are 
discussed  in  reference  (1)  where  a  measure  of  performance  of  mini- 
mization of  maximum  flow  time  is  selected  and  the  problem  solved 
as  a  travelling  salesman  problem.   References  (2),  (3),  and  (4) 
treat  the  related  problem  where  set-up  times  are  sequence-independ- 
ent. 

III.   JOBS  PROCESSED  IN  FCFS  ORDER 

Throughout  this  paper  the  following  notation  is  employed. 

p.  =  known  processing  time  for  job  i  independent  of  sequence, 

i  =  1 , . . . ,n 

d.  =  known  due-date  for  job  i 
l 

s.  .  =  set-up  time  for  job  j  given  that  it  follows  job  i, 
1 » J 

i  =  0,1, ... ,n 
c.  =  completion  time  for  job  i 

It  is  convenient  to  combine  the  terms   p.   and   s.  .   to 

J        i»J 

form  p.  .  ,  the  total  set-up  and  processing  time  for  job   j   given 
that  it  is  immediately  preceeded  by  job  i.   It  is  also  convenient 
to  add  an  extra  job  designated  as  job  zero  to  represent  the  pre- 
liminary idle  condition  of  the  machine.   Thus,   s  .   j  =  l,...,n 
is  the  set-up  time  for  job  j  given  that  job  j  is  processed  first. 


A.   THE  FCFS  ALGORITHM 

The  algorithm  begins  by  ordering  the  jobs  in  the  FCFS  sequence, 
If  there  are  no  late  jobs,  this  sequence  is  optimal.   If  not,  one 
or  more  jobs  must  be  excluded  from  the  sequence  to  produce  an  op- 
timal processing  schedule.   Jobs  once  excluded  from  the  original 
sequence  are  not  considered  again  and  can  be  processed  in  any  order 
after  the  last  job  included  in  the  optimal  schedule.   Suppose  job 
j   is  the  first  late  job  in  the  original  sequence,  then  from  among 
jobs   l,...,j   that  job  is  excluded  which  minimizes  the  time  re- 
quired to  process  the  remaining  jobs.   This  rule  is  repeated  each 
time  an  exclusion  is  required. 

The  steps  in  the  algorithm  are  as  follows: 

1)  Number  the  jobs  according  to  the  order  in  which  they 
must  be  processed  and  form  the   (n+1)  by  n  matrix 

(p.  .).  The  element  in  the  Oth  row  and  the  j th  column 
represents  the  processing  time  plus  the  set-up  time  for 
job   j   given  that  job   j   is  processed  first. 

2)  Augment  the  matrix   (p.  .)   with  an  additional  row  which 
contains  the  due-dates  for  the  jobs. 

3)  Compute  the  completion  time   c,  i=l,2,...k  where 

c.=c.n+p.1.,     c   =  0  . 
l    l-l   rx-l,i        0 

For  the  first  k  such  that   c,  >  d   ,  go  to  step  4.   If 
c.  £  d.  ,  i=l,...,n,   then  the  minimum  number  of  late 
jobs  is  zero. 


4)  For  each  column  l,2,...,k   compute  the  savings   S.   where 

'Si  =  Ci+1  "  ^i-l^i-l,!^ 

5)  From  the  matrix,  delete  row  and  column  h   such  that 

S.  =  max    {S.,  i=l,...,k}  . 
n  1 

Go  to  step  3  and  recompute   c.'s   starting  with  job   h+1. 

(Note  that  after  the  matrix  is  reduced,  the  subscripts 

refer  to  relative  positions  of  columns  in  the  matrix 

rather  than  job  numbers  in  the  computational  formulas 

for   c.   and   S.  .) 
l        l 

The  set  of  columns  in  the  final  matrix  specifies  the  largest 
set  of  jobs  which  can  be  processed  before  their  due  dates.   The  re- 
maining jobs,  which  will  be  late,  can  be  processed  in  any  order 
following  the  last  early  job. 

B.   EXAMPLE 

Let  the  augmented  matrix  of  step  2  be  that  shown  in  figure  1. 


0 
1 
2 
3 

4 


8 


1 

2 

3 

4 

5 

6 

7 

8 

21 

13 

35 

14 

29 

35 

24 

32 

21(20) 

12 

15 

26 

11 

16 

18 

20 

39 

33(23) 

26 

15 

16 

22 

14 

20 

40 

59(57) 

46 

11 

11 

30 

24 

48 

105 

30 

35 

13 

24 

50 

70 

5 
80 

21 
29 

81 

13 
15 

11 

39 

40 

48 

50 

70 

80 

81 

100 

Figure  1:   The  Augmented  Matrix 


The  elements  from  row  9  have  been  moved  up  under  the  main 
diagonal  for  convenience.   Step  3  has  been  carried  out  until 
c„>  d„  .   c.'s   are  shown  in  the  main  diagonal  spaces.   Step  4 
has  been  completed  with   S.'s   shown  in  parenthesis  in  the 
main  diagonal  spaces.   The  quantity   S~   is  max  {S., 1=1,2,3}  , 
so  delete  row  and  column  3  and  return  to  step  3. 
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Figure  2:   The  First  Reduced  Matrix 

The  reduced  matrix  is  shown  in  figure  2  and  again  steps  3 
and  4  have  been  completed  showing  max   {S.,  i=l,2,4,5}   to  be 
S. .   Delete  row  and  column  4  and  return  again  to  step  3.   The  re- 


sulting reduced  matrix  is  shown  in  figure  3. 
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Figure  3:   The  Second  Reduced  Matrix 


Steps  3  and  4  show  now  that   S7   is  maximum,  so  proceeding 


as  before,  delete  row  and  column  7 
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Figure  4:   The  Solution  Matrix 


It  can  be  seen  in  figure  4,  after  deleting  row  and  column  7 


and  carrying  out  step  3  that  now  all  remaining  c,  ^  d.   and  the 


optimal  schedule  is  (1,2,5,6,8).   This  solution  also  has  the  pro- 
perty that  for  all  solutions  in  which  only  3  jobs  are  processed 
late,  this  permutation  completes  processing  of  the  early  jobs  in 
minimum  time. 

C.   PROOF  OF  OPTIMALITY 

As  is  done  in  references  (3)  and  (4),  jobs  will  be  separated 
into  two  disjoint  sets,  E  and  L,  corresponding  to  those  jobs  which 
are  processed  early  and  those  which  are  not  according  to  the  current 
schedule.   Define 

J.   =  i    job  in  FCFS  sequence  of  jobs. 

E.   =  those  jobs  out  of  the  first  k  which  have  been  retained 
in  the  processing  schedule. 

L,   =  those  jobs  out  of  the  first  k  which  have  been  excluded 

k  J 

from  the  processing  schedule. 
|  E,  |  =  the  number  of  jobs  in  set  E,  .   (Similarly  for  |l  |). 

Given  that  out  of  the  first  k  jobs  in  the  sequence   | L,| 
must  be  excluded,  then  it  is  optimal  to  place  those  jobs  in  L, 
which  allow  the   |  E,  |   jobs  in   E,   to  be  completed  in  minimum 
time  since  this  will  allow  the  remaining  n-k  jobs  the  greatest 
opportunity  to  be  completed  before  their  due-dates. 

The  optimality  of  this  criterion  can  be  seen  by  supposing 

that  for  some  current  schedule  with   |L  |   =  j,   2  ^  k  <  n,   the 

next  job  in  the  sequence  has   c.  ,  ,  <  d.  ,  n  .   In  this  case   L.  n  =  j 
J  n  k+1    k+1  '  k+1 ' 

and   c,  , ,   is  minimal.   On  the  other  hand,  if   c.  , .  >  d.  .  ,  then 
k+1  k+1    k+1 


8 


|L    I  =  j+1.   The  job  to  be  placed  in  L   ^   is  job   J   where 

J   is  in  E.  U  (J,  .,}   and  S   =  max   {S.:  J.   in  E,  U  {J,  ,,}}  . 
r  k.  k+1         r  1   1      Tc     k+1 

1 

Then  c.  .._  -  S  £  c.  ■ .  -  S .  ,   J.   in  E.  U  (J,  ,J  •   Hence,   c'  , 
k+1    r    k+1    11       \a  k+1  k 

the  adjusted  time  to  complete  job  k,  given  that  one  more  job  has 
been  placed  in  L,  ,  is  minimized  by  selecting  J  for  placement 
in  L,  .   Optimality  is  thus  proved. 

D.   COMPARISON  TO  MOORE'S  ALGORITHM 

If  s,  .=  s_  .=...=  s   .  =  s.  ,   i.e.,  if  the  set-up 
1,1     2, j  n,j    j 

times  for  all  jobs  are  independent  of  sequence,  then 

p.  .  =  p.  +  s.  .  =  p.  +  s.  ,   is  the  combined  constant 
i>J    J     i»j    J    J 

set-up  and  processing  time  for  job   j .   The  elements  in  the  j 
column  of  the  problem  matrix  will  now  all  be  equal.   Hence, 
S.  =  p.  +  s.   for  each   i   independent  of  sequence.   The  algorithm 
presented  in  this  paper  will  choose  that  job  out  of  the  first   k 
with  the  largest  processing  time  and  place  it  in  L  ,  as  does 
Moore's  algorithm. 


IV.   NO  RESTRICTION  ON  PROCESSING  SEQUENCE 

In  this  formulation  all  n  jobs  are  assumed  to  be  simultan- 
eously available  for  processing  and  no   restrictions  are  imposed 
on  the  processing  sequence.   This  algorithm  uses  a  problem  matrix 
similar  to  that  prescribed  by  step  2  of  the  FCFS  algorithm  just 
described  except  now  both  the  upper  and  lower  triangles  will  be  re- 
quired since  job  numbered   (i+1)   may  preceed  job   i. 


A  branch  and  bound  technique  is  developed  in  which  branching 
to  a  node  labeled   i  means  "process  job   i  next."  The  bound 
computed  at  each  node  is  a  lower  bound  on  the  minimum  number  of 
jobs  which  will  be  processed  late  if  the  solution  set  represented 
by  that  node  is  used. 

A.   THE  BRANCH  AND  BOUND  ALGORITHM 

Steps  1  and  2  are  the  same  as  in  the  FCFS  algorithm  but  are 
repeated  here  for  reference. 

1)  Number  the  jobs  in  non-decreasing  due-date  order  and 
form  the   (n+1)   by  n  matrix   (p..).   The  element  in 
the  0    row  and  the  j    column  represent  the  pro- 
cessing time  plus  the  set-up  time  for  job  j   given 
that  job   j   is  processed  first. 

2)  Augment  the  matrix   (p.  .)   with  an  additional  row  which 
contains  the  due-dates  for  the  jobs. 

3)  At  the  beginning  node  designated  "ALL",  compute  a  first 

bound  B  on  all  solutions  by  comparing  p_  .   to   d.  , 

*J  ♦  3  J 

j=l,...,n.   For  each   i   such  that  p   .  >  d.  ,  increment 

J  Vj     j 

B  by  1  and  delete  the  corresponding  row  and  column. 
(This  step  removes  from  consideration  any  jobs  which 
cannot  be  processed  on  time  no  matter  what  their  position 
in  the  sequence.) 

4)  Branch  from  the  current  node  labelled  k   to  each  remain- 
ing job  as  follows:   For  job   j  ,  subtract   p,  .   from 
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all  elements  in  row   (n+1)   and  compare  the  resulting 
elements  in  row   (n+1),  call  them   (d.(j),   to  the  corre- 
sponding elements  in  row  j.   For  each  p.  .  >  d.(j), 
increment  the  bound  on  node  j   by  1. 
5)   Branch  now  as  described  in  step  4  from  the  node  with  the 
least  bound.   If  there  is  a  tie,  branch  from  the  node  with 

minimum  d..   If  there  still  is  a  tie,  branch  from  the 
l 

node  among  the  tied  nodes  with  minimum  p.  ..   Repeat 
step  5  until  branching  is  completed.   An  optimal  sequence 
is  described  by  tracing  back  through  the  solution  tree. 


B.   EXAMPLE 


The  problem  matrix  is  as  shown  in  figure  5 
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Figure  5:   Problem  Matrix 


Carrying  out  step  3,  obtain  a  first  bound  on  all  solutions 
by  comparing  row  0  to  row  6.   Note  p_   >  d   ,  so  set   B=l   for 

U  «  i        _L 

node  "ALL"  and  delete  row  and  column  1. 
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C.   DISCUSSION 

It  may  appear  that  a  relatively  large  number  of  bounds  must 

be  computed  as  compared  to  well  known  algorithms  (5)  presently 

used  to  solve  the  very  similar  traveling  salesman  problem.   This 

may  turn  out  to  be  the  case.   However,  the  algorithm  takes  great 

advantage  of  one  of  the  trade-offs  available  in  branch  and  bound 

methods  in  that  bounds  are  very  easily  calculated.   The  possible 

requirement  to  compute  more  bounds  then  does  not  necessarily  imply 

reduced  efficiency.   In  fact,  for  the  example  problem  given  only 

"v2     •  2 
2  I      (n-i)   +  n  =  63  additions  and  comparisons  were  required  to 

i=l 
solve  the  problem.   An  n-city  traveling  salesman  problem  solved  by 

the  algorithm  of  Little,  et.  al.,  reference  (5),  would  require 

n-2 
more  than   4   £   (n-i)(n-i-l)  +  3n(n-l)  =  140   additions  and  com- 

i=l 
parisons  even  if  branching  were  only  required  directly  down  a 

single  branch.   The  numbers  of  computations  indicated  above  repre- 
sent lower  bounds  for  each  algorithm.   For  the  algorithm  presented 
here,  if  branching  had  been  required  into  the  second  level  from 
each  of  the  nodes  at  level  1  and  the  remaining  computations  com- 
pleted as  in  the  example,  the  required  number  of  additions  and 

nr2  ,    2 
comparisons  would  have  been  approximately   2(n-l)  I      (n-i)   +  n  =  237 

i=l 
No  firm  conclusion  can  be  drawn  from  these  comparisons,  but  they 

are  presented  as  an  indication  that  the  computational  efficiency 

of  this  algorithm  deserves  further  investigation. 
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